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Path integral and Monte Carlo simulation 
To characterize the thermal state, we can use path integral, 
Z - try End = | DAexp(-5[4)). 
M 
The integration domain M (of an SU(N) gauge theory) is 
M= SU(N) Volume, 


We want a tool to evaluate Z without exponential complexity. 


Monte Carlo method: Consider the case e7314] > 0. 
Generate an ensemble {A;}; following e314, and evaluate 
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More about Monte Carlo simulation 


Do we really circumvent exponential complexity using MC method? 


- Typical values of |O(A;) — (O)| 
JN l 
It indeed solves the exponential complexity for operators satisfying 


(O(A)) 
Typical values of |O(A;) — (O)| 


N 

1 
E f — A; 
rror of = 3 O(A;) 


~ (8+ Volume)~*. 
It has been quite successful to understand Hadron structures, 
thermodynamics of finite-temperature QCD, etc. 


This argument is true only when e S4 > 0. 
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Sign problem and Exponential complexity 


To use Monte Carlo method when e14] 4 0, we generate the 
ensemble {A;}; following the phase-quenched distribution e12e(914). 


(HA) PD He Doe Ole CUD 
(eilm(SIA))) a. Ss e—iIm(S[A]) 


(O(A)) = 


Since (ED ca: _ e78 Volume Af Thus, 
Necessary N of configurations > e2? Volume As, 


Exponential complexity revives due to the sign problem. 


Question 
Can we make Af - 0 by inventing a clever technique? | 
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Method: Path integral on Lefschetz thimbles 
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Lefschetz thimble for Airy integral 
Airy integral is given as 


3 
Ai(a) - [Zoi (= + ar) 


Complexify the integration variable: z = x + iy. 


4 


Integrand on R, and on A 
(a= 1) 
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Multiple integrals on Lefschetz thimbles 
Oscillatory integrals with many variables can be evaluated using the 
“steepest descent” cycles Js: (classical eom S'(z0) = 0) 


f d"™x eS) — wam | dee"), 
M 


o o 


Unlike one-dimensional case, the steepest descent manifold is not 
uniquely defined. 
=> Use of the homology H„ (Mc, {e729 < 1}) becomes quite 


essential: 
> ZIJ), 


XOZ [Ka]. 


K 


H,,(Mc, (ep Ret) < 1} 


Hp (Me (e0 YO < 1}) 


K 


[Pham, 1967; Kaminski, 1994; Howls, 1997] 
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Multiple integrals on Lefschetz thimbles 


Concrete construction Pick up a metric ds? = g,;dz' ® d2/, and 
consider the gradient flow: 


Jo are called Lefschetz thimbles, and Im[S] is constant on it: 


lim z(t) = co} l 
t--00 


Similarly, Ko = {z(0)|z(oo) = ze}. 
[Pham, 1967; Kaminski, 1994; Howls, 1997, Witten, arXiv:1001.2933, 1009.6032] 
[Christoforetti et al. (PRD(2012)), Fujii et al. (JHEP 1310), etc.] 


In = 4200) 
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Monte Carlo algorithm on Lefschetz thimbles 


Monte Carlo algorithms on thimble(s) has been developing: 


e Langevin on J, Cristoforetti, Di Renzo, Scorzato, 1205.3996 


ə Hybrid MC on J> Fujii, Honda, Kato, Kikukawa, Komatsu, 
Sano, 1309.4371 


ə Contraction algorithm Alexandru, Basar, Bedaque, 1510.03258 


ə Generalized thimble Alexandru, Basar, Bedaque, Ridgway, 
Warrington, 1512.08764 


These methods generate ensembles {z;}; on J, following e 


d”zi 
| NO) 
by d"zO(z) exp(-h715(2)) ~ N DDA 1 Jd” z;] i 
YA YA N LY 1 EULA 


N.B. satisfies resurgence for “nice” S (Berry, Howls, '91, Howls '97). 
To 


—S(z). 
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“Reweighting factor” for Lefschetz thimble 


We define the “reweighting factor” of the Lefschetz-thimble approach 
by 


die No Jz, dzexp(—S(z)) 
dig no S7, |dz| exp(—Re(S(z))) 
cf. Reweighting factor in the conventional approach: 


Ju Ax exp(-5(t)) 
Ju Oe exp(—Re(S(x))) 


Idea of this talk: Comparison of these Af tells us the property of the 
Lefschetz-thimble approach. 


exp(--8 : Volume Af) = 


exp(--8 : Volume Af) = 
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Applications 


Applications: Case studies of fermionic sign problem 
e@ One-site Hubbard model (1509.07146, with Tomoya Hayata, Yoshimasa Hidaka) 


ə Multi-flavor massless QED% (1612.06529, with Motoi Tachibana) 
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Applications 


Case 1 One-site Hubbard model 
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AIET One-site Hubbard model 


Path integral for one-site Hubbard model 
We consider the (0 + 1)-dimensional fermion model, 


s= far (0/2 u= iole). 


The path-integral expression for the one-site Hubbard model 


(2 = 3 fy ard(r)): 


p f i 2 B2 
Zan | do (1 + chip t+ ¥/2))? 9897/20, 
2m7U R AB 


Fermion Det 


Integrand has complex phases causing the sign problem. 


p is an auxiliary field for the fermion number density: 
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One-site Hubbard model 


Applications 


Behaviors of number density and Lefschetz thimbles 


Lefschetz thimbles with —0.5U < u < 1.5p: 


EN] 
- 
a 
His 
D 
= 
= 
ey ere ees a: 
Se APTI RRR SR 
Fae: Whe se Re ee Hehe he 
BO BA ee eee ee 
man po NP... ONG Ka z 
a gir YA a a UP 
YA S | 
PAK LEFT H 
AA eee Oe 
TS EE 
asi PAG IA ~ 8 
Ba” Capa ge 
gee a cal en eae de ae ae =e 
a Aa Se ge te i ee ae eel 
ee > ma aaae S- |E 
AA ia manda 5 
os y = mi 2 
aa AI 
x2 a 
WA Pan L fs 
Da Na o a mar wa 
AA AA Ne 
NG a e $C 
ma ape ee ee eS 
DNNN ee ne gir 
SIPON a paagi la 
N BE GA. 7 
NGYN f2 2 ee Se 
TERR or T. 
EEEIEE paa fe 
sl’ 
Aaa Sew ee eee PAA 
2 Mm Pa 3-9 


Number density n with exact result and one-thimble result: 


> 


1.5 JU 
(YT, Hidaka, Hayata, 1509.07146)(cf. Monte Carlo with 1-thimble approx. gives a wrong result: 


1 


0.5 


0 


-0.5 


Fujii, Kamata, Kikukawa,1509.08176, 1509.09141, Alexandru, Basar, Bedaque,1510:03258.) 
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Results 
Results for GU = 30, 100: (YT, Hidaka, Hayata, 1509.07146) 
pU TRD ot i i = - n 
mmal i i if Ts bu =10 PF HA 
--- 3-thimble — exact i a 
— 5-thimble 225 i tible; 2 ol 1 
1 


Í 
1.5) 


1 
\ 


u/U 


i | ni n st $ T3 7 1 1 
- 05 0.5 1.0 1.5 ki 1 1 
$ i 

1 


(cf. Fujii, Kamata, Kikukawa,1509.08176, 1509.09141; Alexandru, Basar, Bedaque,1510.03258) 
Consequence 


In order to describe the step functions, we need interference of 
complex phases among different Lefschetz thimbles. 
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Cit a ales 
Semiclassical partition function 


Using complex classical solutions Zm, let us calculate the semiclassical 
partition function (YT, Hidaka, Hayata, 1509.07146): 


= 1 
Za = > eon = e 50009, (« C KO 5) ad ji 


m=— CO 


This expression is a good approximation for —1/2 < u/U 8 3/2. 
ny Eny 


[U 
-1/2 0 1 Za 
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Sign problem after Lefschetz-thimble deformation 
This computation means that the sign problem exists after the 


Lefschetz-thimble deformation. 
Compute the reweighting factor: 


Ji gy, dee) >” 75m 03 (1679) earl) 
05 ( 


i= 7. no Helga) ~ ST e-Re(Sm) = 0, e-277/BU) 


At u = O, for example, 


Sa In EO ka a 


Te, a) ~ P 


We find nonzero Af. 
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Cire peal 
Comparison between naive reweighting and thimbles 
2.0 


0.0 
-0.5 0.0 0.5 1.0 1.5 


Conseguence 
Lefschetz-thimble method reduces Af in the one-site Hubbard 
model, but it is still nonzero. 
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Applications 


Case 2 Massless QED2 
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Maes GED 
Multi-flavor massless QED, 


2-dimensional U(1) gauge theories with Ny massless fermions: 


Ng 
Z= J DA e54 I DYDyexp | -X ii Pap, [Da - Ha] ba 
a=1 


Since the fermions are massless, the nonzero topological sectors do 
not appear: 

27 o, 27 1 

gt + z inde + *d@ + dr. 


o photon gauge 


A= 


o-dependence is computable using the anomaly equation, and does 
not have the sign problem. 
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HECU Massless QED9 


Toron-field integral of multi-flavor massless QED» 
The toron-field integration becomes (7 = L/: temperature) 


1 27 
Z= i dhodh, exp |--Z Fto, m) A 
0 


where F is the fermion one-loop free energy (pw, = Lua/(27)), 


N 
1 ? T Lo 27 | 1 : 
F=N; (1-5) - EY fin (ear) 


621 n=l 


+ In (2 + gr AN) tn G 4 ee (nth —1+H4)4 a) 
Hin G + er rai) ae \. 


In the limit 7 — O, we can use the mean-field approximation with 
complex saddle points (1612.06529, with Motoi Tachibana) (cf.1504.02979, with 
Hiromichi Nishimura, Kouji Kashiwa). 
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Gradient flow for massless QED» 


All the relevant complex saddle points have real F: 


ZI NAY 
HIN Hs 


aaa 
bone 


YA (Al 
BK NI 
YAI 
zi YA Mi M N NI 

Figure: Gradient flow in the 'Re(hi)- rIm(h 0) plane at Ng -3 
(1612.06529, with Motoi Tachibana) 


a a = = = nag 


Phase structure of multi-flavor massless QED» 


At the zero-temperature with finite L, 7 = 0, the first-order 
transition occurs in this model: 


nı — n2 HaL 
A 27 
6 
4 
2 


a 
1 3 5 

NG no 2 3 3 

(a) Fermion number at N; = 2 (b) Phase boundary at Ny = 3 


(1612.06529, with Motoi Tachibana) (The results are consistent with the exact computation 


given by Lohmayer, Narayanan, 1307.4969) 


a> «a > zr = nag 
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Was CE 
Sign problem after Lefschetz-thimble deformation 


In the Ny-flavor massless QED,, there are only Ny relevant saddle 
points in the limit 7 - 0. 


Since the (complex) mean-field approx. is good in this limit after the 
Lefschetz-thimble deformation, one can show that 


Je. Be Ah e-5(h) 


T 7, hje RT) ~1 


This implies that Af = 0 (at least in the limit 7 - 0). 


Consequence 
Lefschetz-thimble method solves the sign problem of massless QED». | 


Yuya Tanizaki (RIKEN BNL Fermionic sign problem Sep 20, 2017 @ Keio 257 27 
p 


Summary 


ə The sign problem is reviewed from the viewpoint of exponential 
complexity. 


ə Using Cauchy's theorem, one can deform the oscillatory integral 
into the sum of steepest descent integrals on Lefschetz thimbles 


e We consider two examples of the fermionic sign problem. 
In the one-site model, the exponential complexity is not solved, 
but Af is reduced. 
In massless QED», the exponential complexity is eliminated by 
using Lefschetz thimbles. 
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